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We show that the strain-induced attraction between Abrikosov vortices has a non-core contribu- 
tion overlooked up until now. This contribution is an example of the universal mechanism of soliton 
attraction in solids revealed in [Phys. Rev. B 66, 14111 (2002)]. The resulting interaction energy 
is larger than that due to the vortex cores at least by a factor In^ k. Consequently, the non-core 
contribution must be included, for instance, in discussions about the orientation of vortex lattices 
with respect to the crystal axes. It is shown to be also important when interpreting the thermal 
anomalies of the transition between superconducting and mixed states. 
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I. INTRODUCTION 

Attractive interactions between Abrikosov vortices in 
type II superconductors have been a topic of great in- 
terest (see, e.g., Refs. ||, |[ |[ ||]). In particu- 
lar, elasticity-driven attraction between vortices proves 
to have important experimental implications ||, |[ |j. 
However, we shall show that this type of attraction was 
considerably underestimated until now. It is because it 
was assumed that vortex cores, considered as cylindri- 
cal inclusions of the normal phase, were the only sources 
of strain. In Ref. Q a universal mechanism of strain- 
induced attraction between solitons has been revealed. 
This is operative even neglecting the strain induced by 
the soliton (vortex) cores (see below). In this Letter 
we show that this non-core mechanism gives the main 
contribution to the elasticity-driven attraction between 
Abrikosov vortices if the Ginzburg-Landau parameter k 
is large. 

We will consider elastically isotropic media only. Al- 
though this case is of some interest per se, its main 
convenience is that it represents a fairly simple case in 
which the leading role of the non-core contribution to 
the elasticity-driven attraction can be shown. With this 
result in mind it is natural to expect that in elastically 
anisotropic media, which have more practical interest, 
the same situation takes place (see Appendix). 

We will discuss several consequences of this mechanism 
of attraction. First of all we will show that its contribu- 
tion to the energy of the vortex lattice (VL) differs from 
that associated with the cores by a factor In^ k or even 
more. Previous authors |[ ^ considered only the core 
contribution so they strongly underestimated this part of 
the VL energy of high-K superconductors. Despite this 
underestimation, the observed correlations between the 
VL and crystal lattice has been successfully explained by 
virtue of the strength of the strain-induced interaction in 
the considered materials. However, for other materials 
which can be studied in future it could be necessary to 
take into account both core and non-core contributions. 



The existence of any attraction between vortices im- 
plies the discontinuity of the transition between super- 
conducting and mixed states. As we will show, the dis- 
continuity of the flux density associated with the strain- 
induced attraction between vortices is strong enough to 
be measurable in typical superconductors. There exist as 
well a latent heat associated with this discontinuity. It is 
expected to reveal itself experimentally as an anomaly in 
the specific heat. This anomaly should be taken into ac- 
count when interpreting the experiments on the specific 
heat of mixed-state superconductors. 

The nature of strain-induced interactions between vor- 
tices can be easily illustrated considering a system with 
an infinite shear modulus (/i — oo) 0. The only pos- 
sible deformation of such a medium is a homogeneous 
dilatation u. The finite size of the sample is, in fact, 
taken into account by introducing this homogeneous de- 
formation and considering a free crystal, i.e. introduc- 
ing implicitly boundary conditions Suppose that a 
density of vortices n is created. The vortex self-energy 
depends, naturally, on u. Taking the state without vor- 
tices as the non-deformed one and omitting for a while 
the repulsive interaction between them, we present the 
change in the energy of the system per unit volume as 
F{u) ~ n{Eo + Eiu) + where K is the bulk mod- 

ulus and the terms in parentheses represent the above 
mentioned vortex self-energy. Minimizing with respect to 
u one obtains the equilibrium deformation of the sample: 
Uoq = —nEi/K. Then, the change in energy becomes 
F(uoq) = nEo — n'^El/{2K), where the second term rep- 
resents the strain-induced attraction between vortices. 
Evidently one cannot forget that some vortex repulsion 
exists as well, a repulsion that provides a finite value of 
the equilibrium vortex density. Below we reproduce this 
result in more detail. 

II. STRAIN-INDUCED ATTRACTION 

We will start by showing that even neglecting the vor- 
tex cores, i.e. neglecting all sources of strain previously 
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considered in Refs. H, |[ there is a strain-induced 
attraction between vortices. Within Ginzburg-Landau 
theory the free energy density can be written as (see, 
e.g., Ref. @) 
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Here ^' is the order parameter, a = a{T — Tc) \s the only 
temperature dependent coefRcient where Tc is the critical 
temperature, and v is the volume of the system. The elas- 
tic degrees of freedom are taken into account in Eq. ( pb| ) 
where Uik is the strain tensor and summation over dou- 
ble indices is implied (see, e.g., Ref. [0). The first term 
in Eq. (|2b|) is responsible for the pressure dependence of 
the critical temperature and the change of the bulk mod- 
ulus corresponding to the normal-superconductor phase 
transition. 

The simplest way to work out the strain-induced vor- 
tex attraction is following the method given in Ref. 
To minimize the free energy over the elastic degrees of 
freedom we distinguish between homogeneous and inho- 
mogeneous deformations [pl|: 



Here, uf^ represents the tensor of homogeneous defor- 
mations and Mi(k) the components of the displacement 
vector in Fourier space. Minimization of Eq. ( [2b| ) with 
respect to all elastic degrees of freedom gives 
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erage. The contribution of the first term reduces to a 
renormalization of the coefficient b in Eq. (pa[). This 
renormalization disappears in the limit fi oo. Due 
to the second term the free energy becomes a non-local 
functional. Note that this non-locality remains as long 
as the shear modulus does not vanish. 

Further minimization of Eq. (^) with respect to ^' 
is not straightforward due to its non-locality. However, 
if = CX3 there is another minimization procedure that 
avoids the treatment of non-local equations. We first 
consider this case and after that we return to fi ^ oo. 



As we have mentioned, in this case the only possible 
deformation is a homogeneous dilatation u. Therefore, 
the free energy (|l|) can be written as 
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where a[u) = a + ru and u is a variational parame- 
ter. Fixing for a while this parameter, i.e. consider- 
ing momentaneously a clamped sample, the form of the 
equations of equilibrium reduces to that of the Ginzburg- 
Landau equations Solving them one obtains, in par- 
ticular, the free energy density close to the transition 
between superconducting and mixed states in terms of 
the magnetic induction B = AttB. For triangular VL's in 
high-K superconductors (Iuk ^ 1) it can be written as 



F=Fs{u) + ^u^ + B[H,i{u) 
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where v — 3^/^(7r/2)^/^. The first two terms, being 
Fs{u) = —a^{u)/ {2b), represent the free energy den- 
sity in the superconducting state at 7J = 0. The third 
term is proportional to the vortex self-energy, where 
Hci{u) is the magnetic field at which this self-energy 
changes its sign. It is known that for high-K supercon- 
ductors the vortex self-energy comes mainly from contri- 
butions of the non-core region So taking Hci{u) = 
{(/)o/[47rA^(it)]} In K, where 4>o is the flux quantum and 
A(m) = {mc^6/[87re^|a(M)|]}^/^ is the penetration length 
of the magnetic field, we shall reveal effects associated 
with these non-core contributions. The last term repre- 
sents the repulsive interaction between vortices that takes 
place at low flux densities. Here Bo{u) — (/)o/[27r\/3A^(u)] 
defines the reference flux density. 

Let us now consider free samples. Then Eq. (^ has 
to be minimized with respect to u. The equilibrium de- 
formation in the superconducting state is Us = ar/{bK), 
where K = K — r^/b. In the mixed state there is, in addi- 
tion, a deformation Um — u^Ug as a result of the creation 
of vortices. Since it is small close to the transition be- 
tween the superconducting and the mixed states, only 
lowest order terms are relevant and the Wm-dependence 
of the repulsion term can be neglected in Eq. (|^). Thus, 
minimizing Eq. (^ with respect to Um we obtain 
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where 7 — K/K and S = Tx\r^ j {Kb)\K ^ In^ k for high-K 
superconductors, i.e. In k = 2k Hci/Hc2 H (here and 
in what follows the values i/ci, Hc2: etc. are refereed 
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to the non-deformed state, u = 0, if it is not exphc- 
itly indicated). The last term of Eq. (|^) represents the 
strain-induced attraction between vortices. Mention that 
it disappears if the shear modulus goes to zero, i.e. it 
is associated with the solid-state elasticity. Let us recall 
that evaluating the vortex self-energy the vortex core was 
neglected. 

The free energy could be obtained, in principle, 
from Eq. with its coefficients corresponding to /x = 
oo, i.e. rV(2K4/3) = and 4^/(37^4/3) = 1- Note 
that there is no essential difference between the func- 
tional form of Eq. (^) for infinite and finite ^. So we 
conclude that the free energy density of aiw isotropic 
type II superconductor has the form of Eq. (0) with the 
renormalized constants h' ~ h — /K^^iy, and {r^ /K)' ~ 
{ryK)[A^,/{iK,„)]. 

III. CORE CONTRIBUTION 

We will deduce the core contribution to the strain- 
induced attraction for ^ = oo in order to compare it with 
the non-core one Following Ref. we model the 
vortices as cylinders of radius ^ (coherence length) of nor- 
mal phase inside a superconducting medium. Let us first 
consider a clamped superconducting medium. To accom- 
modate a normal cylinder inside this medium, the cylin- 
der should be deformed because of the difference between 
specific volumes of normal and superconducting phases 
Vn,s- Such a deformation is simply uq = {Vn — Vs)/Vs 
(if ^ = oo only homogeneous deformations are possible). 
Let us now consider a free sample designating as n the 
density of cyHnders (vortices). The elastic part of the 
free energy density can be written as 



F~mrf — {u-uof + —u^, 



(8) 



where u is the deformation of the sample as a whole and it 
has been taken into account that the bulk moduli of both 
normal and superconducting phases are approximately 
equal (K). Minimizing Eq. (^ with respect to u we 
obtain the equilibrium deformation of the sample: Um — 
nTT^^ug. Therefore, the equilibrium free energy results to 
be 
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The second term of Eq. (H) represents the attraction 
between vortices due to the core-induced strain. Taking 
into account that n — B/(j)Q and (see, e.g., Ref. [Q) 
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where = A-ko? /h and (dTc/du) — r/a, this second 
term can be written as 
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For high-K superconductors one has 2k^ Hd/ Hc2 = Iuk. 
Thus, the ratio between the core contribution to the 
strain-induced interaction and the non-core one [see Eq. 
(0)] is 
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This ratio can be calculated more consistently close 
to Tc- In this region vortex cores can be success- 
fully described within the Ginzburg-Landau theory, i.e. 
no model is necessary here. One has Hci{u) — 
{0o/[47rA2(u)]}(lnK 0.08) in Eq. (§) (see, e.g., Ref. 
1^). Further minimization of Eq. (H) gives a free energy 
of the form Eq. (|^) where the resulting coefficient 5 in- 
cludes both core and non-core contributions. It is such 
that 
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As we see, the model of Ref. Q strongly overestimate 
core contribution in this region. In any case, the strain- 
induced attraction energy due to the vortex cores results 
to be at least one order of magnitude lower than the non- 
core one for high-K superconductors. 

The latter also shows that the contribution to the at- 
traction energy due to the core-induced strain can be un- 
derstood as a correction of the non-core one. Both con- 
tributions depends on orientation of the VL with respect 
to the crystal axes if one takes into account some elas- 
tic anisotropy, being this dependence the same (see Ap- 
pendix). The ratio between both contributions remains 
given by Eq. (^. In Refs. |||, [H § it was suggested that 
the above mentioned dependence of the strain-induced in- 
teraction can explain the correlations between the vortex 
and crystal lattices. Mention that in Rcfs. |j this was 
asserted carrying out a detailed comparison between the 
London and the core- induced elastic energies of the VL's 
in NbSe2. The experimentally observed VL in NbSe2 
does not correspond to the lowest energy one if only the 
London energy is taken into account. But including both 
London and elastic contributions, the experimentally ob- 
served VL coincides with the lowest energy one. This 
coincidence is a luck because the most important contri- 
bution to the elastic energy of the VL's, i.e. the non-core 
one, was overlooked. It might not be the case in other 
compounds ||l^ , in which the non-core contribution could 
be essential. 

IV. THERMAL ANOMALIES 

As a result of the attraction between vortices, the tran- 
sition from the superconducting to the mixed state has 
a discontinuous character. The conditions of minimum 
and continuity of the free energy determines the jump 
of the magnetic induction and the transition magnetic 
field. The attraction does not affect substantially the 
transition magnetic field which remains to be close to 
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Hci- The jump of the magnetic induction AB is such 
that 

AB/Bo ~ hi"2 i^[r^/(^Kb)]K'^ In^ k| , (14) 

where Bq — AttBq. The ratio r'^/{Kh) is the rel- 
ative change in the bulk modulus at the normal- 
superconductor phase transition. Typically its order of 
magnitude is r^/{Kb) = lO-^-^lQ-^ § |l|. Thus, taking 
into account that k-Muk ~ 10"^-;- 10"^ ^ 10- 10^), 
the jump in the magnetic induction is expected to be 
AB ~ 10~^ 10'^ Bq. These values are high enough to 
be experimentally appreciable. 

There is a latent heat associated with this discon- 
tinuous transition. From Eq. (0) it can be written 
as Q ~ TAB{dHci/dT). The contribution due to 
the strain-induced attraction can be roughly estimated 
as ~ 1 mJ mol~^ in NbSe2 and ~ 0.1 mJ mol~^ in 
YBa2Cu307_5. Since this phase transition is normally 
smeared, this latent heat will manifest itself as a specific 
anomaly. Having in mind the order of magnitude of the 
observed specific heat (~ 10 mJ mol~^ in NbSe2 

and ~ 1 mJ mol~^ in YBa2Cu307_5, see e.g. Refs. 
[ p^ , |l5j) one might conclude that this latent heat should 
be observed. 

Let us mention that measurements of the field- 
dependent specific heat have been proposed as a way 
to distinguish between s-wave and d-wave superconduc- 
tors. The specific heat has, in both cases, contributions of 
terms which are proportional to the quasiparticles den- 
sity of states localized in the vortex cores. Assuming 
that the vortex density depends linearly on the magnetic 
field, these terms are expected to be oc 7J in s-wave su- 
perconductors |l^ . In d-wave ones a weaker cx y/H field 
dependence is expected as a result of the quasiparticle 
delocalization [0. In consequence, s and d-wave super- 
conductors are hopped to be distinguishable by presence 
or absence of curvature in the specific heat dependence 
on the external field. 

In Refs. 0, 

m it was pointed out that, even in s-wave 
superconductors, interactions between vortices could in- 
duce this curvature: these interactions provoke a non- 
linear dependence of the vortex density on the external 
field. Only the repulsive interaction was explicitly men- 
tioned, but our results indicate that the attractive inter- 
action could be important as well. Taking all interactions 
into account, the specific heat anomaly associated with 
the transition between the superconducting and mixed 
states can be obtained from Eq. (^ in a similar way to 
that shown in Ref. p9[ |. 

Indeed, the dependence of the specific heat on the mag- 
netic field might be affected by the smeared latent heat 
of the superconducting-mixed state transition. This con- 
tribution can not be quantified exactly, but the orders 
of magnitude obtained above indicate that it could be 
a significant contribution. As we see, when distinguish- 



ing between s-wave and d-wave superconductivity by way 
of measurements of specific heat anomalies much care 
should be taken. 

V. CONCLUSIONS 

We have shown the leading role of the non-core 
contribution to the strain-induced attraction between 
Abrikosov vortices in high-K superconductors. In the 
elastically isotropic case studied, the strain-induced in- 
teraction is due to finite size effects what previous authors 
neglected. The importance of the non-core contribution 
is also expected in a general case. Thermal anomalies 
associated with the transition between superconducting 
and mixed states have been discussed. The resulting 
jump in the magnetic induction and the latent heat of 
this transition was found strong enough to be measur- 
able. The later might be important when interpreting 
experimental data on the specific heat of mixed-state su- 
perconductors. 
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APPENDIX 

A full calculation of the strain-induced interaction be- 
tween vortices in elastically anisotropic media is quite 
complicated and, to the best of our knowledge, has never 
been performed. The usual method to estimate this inter- 
action consists in: (i) to calculate the interaction between 
two vortices in an infinite medium and (ii) to sum it over 
all vortex pairs (see, e.g., Ref. Q). In what follows we 
shall analyze the first step of such a estimations. 

For elastically anisotropic crystals, instead of Eq. ( pb| ) 
one has 

Fel = ^ J ^rifcMifc|*P + ^XiklmUtkUhn^ dv . (15) 

The minimization with respect to the elastic degrees of 
freedom can be worked out as before (/i 7^ 00). 

Let us consider two well separated vortices, one at p = 
and another at p = Pq {po ^ A), directed along the 
z-axis. The strain-induced interaction between them can 
be written as [Q 

F!riPo)^- J -^e'^-for,^r,f,G,,k^kp\fik)\^ (16) 

where — Xaiii^pkakp (Greek subscripts acquire only 
X and y values), and /(k) is the Fourier transform of the 
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difference between the value of the order parameter far 
from the vortices and the square of the order parameter 
amphtude associated with one vortex. 

In Refs. Q it was only considered the vortex core 

using the above described model of normal cylinders: 



(17) 



where = \a\/b. Let us neglect this contribution, fo- 
cusing our attention in the non-core one. Then, because 
the order parameter amplitude is such that H 
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Due to the separation between vortices, the most impor- 
tant contribution to Eq. ( |l^ ) comes from k ^ A~^. For 
these values, the argument of the Bessel function in Eq. 
( |l9| ) is fcp <C 1 inside of the integration interval, so 



non^core 
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dp ^ 271^^-^1111 K. (20) 



As we see, the non-core contribution to -Ff""^ prevails 
over the previously reported core one by virtue of the 
high value of the Ginzburg-Landau parameter k. 

Mention that close to Tj, the expression (p^ ) can also 
be used in order to estimate the core contribution, i.e. 
no model is necessary here. Thus one obtains 



j core (k) 
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(21) 



Therefore, the ratio between core and non-core contribu- 

el 



tions to iT*""" results to be ~ 1/(2 In^ k) (the model of 



normal cylinders slightly overestimates the core contri- 
bution in this region). 
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